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Abstract 

We consider n one-dimensional Brownian motions, such that n/2 
Brownian motions start at time t = in the starting point a and end at 
time t = 1 in the endpoint b and the other n/2 Brownian motions start 
at time t = at the point —a and end at time t = 1 in the point —6, 
conditioned that the n Brownian paths do not intersect in the whole 
time interval (0, 1). The correlation functions of the positions of the 
non-intersecting Brownian motions have a determinantal form with a 
kernel that is expressed in terms of multiple Hermite polynomials of 
mixed type. We analyze this kernel in the large n limit for the case 
ab < 1/2. We find that the limiting mean density of the positions of 
the Brownian motions is supported on one or two intervals and that 
the correlation kernel has the usual scaling limits from random matrix 
theory, namely the sine kernel in the bulk and the Airy kernel near the 
edges. 
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Figure 1: Non-intersecting Brownian motions when ab > 1/2. Here we have 
chosen a = 1 and b = 0.7. For ab > 1/2 we have that the two groups of 
Brownian motions remain separated during the full time interval (0, 1). 



1 Introduction 

Consider n one-dimensional Brownian motions conditioned not to intersect 
in the time interval (0, 1). We assume that n is even and that n/2 Brownian 
motions start at the position a > at time t = and end at b > at time 
t = 1, while n/2 Brownian motions start at —a and end at —b. If we let 
n — > oo and at the same time rescale the variance of the Brownian motion 
with a factor 1/n, then the Brownian motions fill out a region in the tx-plane 
that looks like one of the regions shown in Figures [TH3J 

The shape of the region depends on the product ab. For ab greater than 
some critical value (which in the units that we will be using is 1/2) the 
starting positions ±a and the end positions ±6 are sufficient apart so that 
the two groups of paths remain separated. The Brownian motions then fill 
out two ellipses as can be seen in Figure [TJ Each group essentially behaves 
like n/2 non-intersecting Brownian motions with a single starting and end 
position. This is a variation of Dyson's Brownian motion for the behavior of 
the eigenvalues of a Hermitian matrix whose elements evolve according to a 
Brownian motion [18J. In that model, it holds that at each time t S (0, 1) the 
positions of the paths are distributed like the eigenvalues of a GUE matrix 
and so follow Wigner's semi-circle law as n — > oo, see |27j . For the present 
model and ab > 1/2, we then also expect to find at each time t 6 (0,1) 
the scaling limits for the correlation functions that are known from random 
matrix theory, and that are expressed in terms of the sine kernel in the bulk, 
and the Airy kernel at the edge points. 

The behavior is different for ab less than the critical value 1/2. In that 
case the starting and end positions are not that far apart and the two groups 
of paths will interact with each other. There are two critical times t C) \ and 
i Cj 2 such that the two groups are separated up to the first critical time t C] \. 
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Figure 2: Non- intersecting Brownian motions when ab < 1/2. Here we have 
chosen a = 0.4 and b = 0.3. For ab < 1/2 we have that the two groups of 
Brownian paths starting at ±a come together at the first critical time £ C) i, 
they merge and continue as one group until the second critical time t C: 2, 
after which they split again and end at ±6. 



At i Cj i they merge and continue as one group of paths until the second 
critical time t c 2 when they split again, see Figure [2j This is an extension of 
the case of two starting positions and one end position which was studied in 
[U [7] with the use of Riemann-Hilbert techniques and in [29] by classical 
steepest descent techniques. There it was found that for each time t € (0, 1), 
the scaling limits are still expressed in terms of the sine kernel in the bulk 
and in terms of the Airy kernel at the edge, except at the cusp point at the 
critical time. At the cusp the scaling limits are expressed in terms of Pearcey 
kernels that were first described by Brezin and Hikami [TU] . Tracy 
and Widom [29] also identified a Pearcey process that is further discussed 
in [U [H [25] as well. 

Returning to the model with two starting positions and two end positions 
we have a critical separation ab = 1/2. Then the starting and end positions 
have a critical separation so that the two groups of paths just touch in 
one point, see Figure [3l Here we expect new critical behavior that can be 
expressed in terms of an as yet unknown kernel. 

It is the aim of the present paper to treat the case ab < 1/2, t ^ i c lj *c,2 
with the methods of [HE1IZ]- That is, we use a steepest descent analysis of 
a relevant Riemann-Hilbert problem that was given by the first two authors 
in [p3]. The Riemann-Hilbert problem has size 4x4 and its solution is 
constructed out of multiple Hermite polynomials of mixed type. 

As a result of the asymptotic analysis we find the usual sine kernel in the 
bulk and the Airy kernel at the edge points, see Theorems 12 . 21 and 12.31 below, 
thereby providing further evidence for the universality of these kernels. We 
are not aware of a double integral representation of the finite n correlation 
kernels, so that an asymptotic analysis of integrals as in [29] may not be 
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Figure 3: Non-intersecting Brownian motions when ab = 1/2. Here we have 
chosen a = 1 and b = 1/2. For ab = 1/2, we have that the two groups of 
Brownian paths touch each other at a critical time t c . 



possible in this case. 

The connections between random matrices and non-intersecting Brow- 
nian motions are well-known, see for instance [20] and [30] for a very re- 
cent paper on non-intersecting Brownian excursions. Unfortunately we do 
not know if there exists a corresponding random matrix model for the case 
studied in this paper. 

2 Statement of results 
2.1 Correlation kernel 

It follows from the classical paper of Karlin and McGregor [21| that the 
(random) positions of the Brownian motions at time t S (0, 1) are a deter- 
minantal point process. This means that there is a kernel K n so that for 
each m we have that the m-point correlation function 

n! f f 
(n-my.J J 

where p Ht t(xi, . . . , x n ) denotes the joint probability density function for the 
positions of the paths at time i, is given by the determinant 

Rm (•£].>•••) %m) — &£fci(K n {xi , Xj ))i I j'=l ) ... J m) 

see [25] . 

Indeed, we have by [21 J that 

p n> t(x 1} ...,x n )(x det(P n (i,ay,Xjfc)) J> fe = i i ... )n det(P n (l-i,x fe ,6j)) Jifc= i > ... )n 

(2.1) 
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if the non-intersecting paths start at a\ < a>2 < ■ ■ ■ < a n at time t = and 
end at 61 < 62 < • • • < 6 n at time t = 1, where P n (t,a,x) is the transition 
probability for Brownian motion with variance 1/re, 



Thus (12. ip is a biorthogonal ensemble [8], which is a special case of a determi- 
nantal point process. These ensembles are also further studied in connection 
with random matrix theory in [T7]. The correlation kernel is obtained by 
biorthogonalizing the two sets of functions P n (t,aj,-) and P n {l — t,-,bj) 
which results in functions <f>j and ipj for j = 1, . . . , n, say, and then putting 



n 



K n {x,y) = ^<t> j {x)^{y). (2.2) 

3=1 

In the situation of the present paper we have to take the confluent limit 
dj — > a, bj — ► b for j = 1, ...,n/2, and aj — > —a, 6j — > —6 for j = 
n/2 + 1, . . . , n. Then we continue to have a biorthogonal ensemble and the 
structure of the kernel (|2.2p remains the same. In [13] the functions <f>j and 

are written in terms of certain polynomials that were called multiple 
Hermite polynomials of mixed type. We will use here only the fact that K n 
is expressed in terms of the solution of a Riemann-Hilbert problem, see |13] 
and (f2~TB|) below. 

We assume throughout that 

0<ab<l/2 and 1 6 (0, 1) \ {t Cjl , i C)2 }. (2.3) 
So we exclude the critical times from our considerations. 

2.2 Limiting mean density 

Our first result deals with the limiting mean density 

p(x) = lim -KJx,x). (2.4) 

n^oo n 



Theorem 2.1 Let a, b > and t G (0,1) such that (I2.3p holds. Then the 
limiting mean density (|2,4p of the positions of the Brownian paths at time t 
exists. It is supported on one interval [—2:1,2:1] if t £ (t c ,i,tc,2) and on two 
intervals [—21, —22] U [22, 21] if t G (0, t Cj i) U (i C) 2, 1)- In all cases the density 
p(x) is expressed as 

p(x) = -\lmC(x)\, (2.5) 
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where £ = £(x) is a solution of the equation 



. i 2: ; . / rr a 2 Ir 1 2 



^ i(l-i) e + U (1 _ i)2 ^ (i_ t )a + t (i_ t )J^ 

/ 2fr 2 1 \ 6 2 z 2 1 J 

+ Vt(l-t) 3 t 2 (l-t) 2 )^ t 2 (l-t) 4 

TTze density p is real and analytic on the interior of its support and it van- 
ishes like a square root at the edges of its support, i.e., there exists a constant 
c\ such that 



p(x) = — \x =F (1 + o(l)) as x — > ±zi, x £ suppp, (2-7) 

7T 



c l I -,- 1 1/2 

— \x =F 

7T 

and, in case i G (0, i Cj i) U (i C) 2, 1), there exists a constant c 2 such that 

p(x) = — \x =F Z2\ 1 ^ 2 (l + o(l)) as s -> ±Z2, s S suppp. (2. 
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2.3 Scaling limits of the kernel 

As in [1] and [6], the local eigenvalue results are formulated in terms of a 
rescaled version of the kernel K n , 

K n (x,y)=e n ^- h ^K n (x,y) (2.9) 

for some function h. Note that this change in the kernel does not affect the 
determinants 

det (K n ( 

that give the correlation functions for the determinantal process. 

Theorems 12.21 and 12.31 show that the kernel has the scaling limits that 
are universal for unitary ensembles in random matrix theory [15], namely 
the sine kernel in the bulk, and the Airy kernel at the edge. 

Theorem 2.2 Let a, b > and t £ (0, 1) such that (12. 3p holds. Let Z\, z 2 be 
as in Theorem 12. 1 1 Then there exists a function h such that the following 
holds for the rescaled kernel (12.9ft . For every xo £ (— z±,Zi) (in case t £ 
(tc,l,t c ,2)) or xq £ (-2i, -z 2 ) U (z 2 , Z\) (in case t £ (0, i C) i) U (i C)2 , 1)) and 
ii,ti£l, we have 

1 - / u v \ smir(u-v) 
lim ——K n x + —f-^r, x + -— ^ = ^ -A 2.10 

- np(x ) V np{x Q ) np{x ) ) ir(u - v) 
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Theorem 2.3 Let a, b > and t € (0, 1) such that (12. 3p holds. Let z±, Z2 be 
as in Theorem 12, 11 Then there exists a function h such that the following 
holds for the rescaled kernel (12.91) . For every we have 



n-^oo (cin) 2 / 3 



u 



,Zi + 



(cin) 2 / 3 ' (cin) 2 / 3 
Ai(«) Ai» - Ai'(u) Ai(v) 



u — v 



(2.11) 



where Ai is the usual Airy function, and c\ is the constant defined in ([27 
In addition, if t £ (0,i C) i) U (t c ,2> 1); then for every u,v £ R, we have 



lim 



1 



u 



,Z 2 



n-*oo (c2?l) 2 / 3 \ (C2n) 2 / 3 ' (C2n) 2 / 3 

Ai(u) Ai'(u) - Ai'(n) Ai(v) 



u — V 



(2.12) 



where c 2 is the constant appearing in (12. 8|) . 
Similar results hold near —z\ and —z 2 - 

2.4 The Riemann-Hilbert problem 

Our proofs are based on the fact that the kernel K n can be written in terms 
of the solution of the following Riemann-Hilbert (RH) problem, see |13j . We 
look for a matrix valued function l":C\R^C 4x4 such that 

RH problem for Y 

(1) Y is analytic on C \ R, 

(2) for x € R, it holds that 



Y + (x) = Y_(x) 



(I Wl >1 {x)W2,l{x) W 1)1 (x)w 2; 2(x)\ 

1 wi >2 (x)w2 J i(x) wi )2 (x)w 2;2 (x 

1 

\0 1 / 



(2.13) 



where Y + {x) (Y_(x)) denotes the limiting value of Y(z) as z approaches 
x from the upper (lower) half plane, 

(3) as z — > oo, we have that 

/z n / 2 \ 

z n / 2 
Y(z) = (I + 0(l/z)) " \ " " . (2.14) 



, z nl1 








^ 





z n/2 














z -n/2 





V o 








z -n/2j 
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In (|2.13p we used the functions 



W2,l(x) 



-f t (x 2 -2ax) 



w 1)2 (x) 



e -f t (x 2 +2ax) 



s ^ iy (x a +2te) 



(2.15) 



This RH problem can be seen as a generalization of the RH problem for 
multiple orthogonal polynomials |31| . which in turn is a generalization of 
the RH problem for orthogonal polynomials [19J. 

The kernel K n then takes the following form, see [13] : 



K n (x,y) 



1 



2m{x — y) 



(o o io 2il (y) w2,2(y))Y+ 1 (y)Y+(.x) 



w 1)2 (x) 



(2.16) 



/ 



As explained in [13] the solution of the above RH problem is unique and is 
built out of polynomials that satisfy certain orthogonality conditions that are 
a combination of the conditions satisfied by multiple Hermite polynomials 
of type I and II. Therefore they were called multiple Hermite polynomials of 
mixed type. The multiple orthogonal polynomials of mixed type and their 
RH problem are related to multi-component KP as shown in [3]. 

The formula (|2.16p may be seen as a Christoffel-Darboux formula since 
it equates the sum (|2.2p with the right-hand side of (|2.16p which is of the 
form 



fi(x)gi(y) + f2{x)g 2 {y) + h{x)gz{y) + Mx)gj{y) 

X — V 



x-y 

with only four terms fj(x)gj(y) in the numerator. A similar Christoffel- 
Darboux formula was derived in [T2] for the case of multiple orthogonal 
polynomials. 

We use the Deift/Zhou steepest descent method for Riemann-Hilbert 
problems and apply it to the Riemann-Hilbert problem stated above. This 
yields strong and uniform asymptotics for Y as n — > oo and then also for 
K n due to the formula (I2.16|) . 



2.5 Riemann surface 

As in [3l[6ll7l[23|l26j the asymptotic analysis of the Riemann-Hilbert problem 
is based on a suitable Riemann surface that in our case is given by the 
equation (|2.6|) . The determination of the equation for this surface is more 
complicated in this case, since we could not find an explicit differential 
equation that is satisfied by the multiple Hermite polynomials of mixed 
type. We found the equation (|2.6p only after numerical experimentation 
with Maple. 
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The equation (12, 6\\ has 6 branch points. For ab < 1/2 there are two 
critical times < t c \ < t C) 2 < 1 such that for t Cj % < t < t c ^ four branch 
points are purely imaginary and two branch points are real, while for < 
t < t c i and t C 2 < t < 1, four branch points are real and two branch points 
are purely imaginary. We prove this in Section 13.31 

The rest of the paper is organized as follows. In Section [3] we give more 
details about the Riemann surface associated to (|2.6|) that will be used in the 
asymptotic analysis of the Riemann-Hilbert problem for Y, defined in (|2.13p 
and (|2.14p . In Section 0] we use the Deift/Zhou steepest descent method to 
analyze this Riemann-Hilbert problem. It leads to the proofs of Theorems 
l2"3HO in Section [5j 

3 The Riemann surface 

In this section, we study the Riemann surface associated with the equation 
(|2.6p which will be used in Section 0] for the asymptotic analysis of the 
Riemann-Hilbert problem. When we talk about the Riemann surface we 
will always mean the compact surface that arises after resolution of the 
singularities of (|2.6[) . 

3.1 Nodal singularities and genus 

First of all we note that (|2.6p is singular, since clearly both partial derivatives 
vanish at £ = z = and so the origin is a singular point. There are also 
singular points at infinity and to study those it is convenient to introduce 
homogeneous coordinates [z : £ : w] and the homogeneous equation 

,4 2 n 1 22 (a 2 b 2 1 \ a2 2 

{2b 2 1 \ . 2 b 2 2 2 n 

(3.1) 

in projective space P 2 . 

Proposition 3.1 The equation (13. ip has three nodal singularities at [0 : : 
w], [z : : 0], and [z : t (^_ t ^ ■ 0]. Each of these singularities corresponds to 
two points on the Riemann surface, and the Riemann surface has genus 0. 

Proof. To verify that [z : : 0] is a node, we set z = 1 in (|3.ip and we 
note that the resulting equation in £ and w has vanishing partial derivatives 
at £ = w = 0. The quadratic part (in £ and w) is 

1 t 2 _ b 2 2 

t 2 (i-t) 2 ^ t 2 (i-ty w 
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which gives rise to two distinct tangents 

£ = ±j^-™ (3-2) 

and therefore [z : : 0] is a node. Similarly, [0 : : w] is a node. 

The third singularity at [z : j^b^y : 0] can De readily seen if we rewrite 
flHZED as 



t(i-t)/ t 2S t(i-ty V 



6 2 



(1 _ t) 2 t(1 _ t) 



(3.3) 



Then replacing z by it = £ — a , setting £ = 1 and taking the quadratic 



part (in w and u) we get u 2 — ^rw 2 which again gives two distinct tangents 
u = ±| w. Therefore [z : trnbj : 0] is a node as well. 

After resolution each of the three nodes leads to two points on the Rie- 
mann surface. 

The genus formula (Plucker's formula |24} Proposition 2.15]) says that 

g = {d- l)(d-2)/2 - k 

for the genus g of a surface of degree d with nodes and no other singulari- 
ties. In the present case we have d = 4 and k > 3. Since (d— l){d — 2)/2 — k 
is always an upper bound for the genus, we conclude that the Riemann sur- 
face has genus zero. We also conclude that there are no other singularities 
besides the three nodal singularities we already found. □ 



3.2 Rational parametrization 

We find a rational parametrization of the Riemann surface by intersecting 
the conic 

f ~ t r 1 _ t f Z + Pt w + Q ZW = ( 3 - 4 ) 

with the equation (13. ip . By Bezout's theorem there are 8 intersection points 
in P 2 if we count intersections according to their multiplicities. It is easy to 
see that the conic intersects equation (|3,ip at the three nodes [0 : : w] , [z : 
0:0], and [z : t ^_ t ^ ■ 0] (see Proposition I3.1|) for any choice of parameters 
p and q. This accounts for at least 6 intersection points, since each of the 
nodes counts at least twice. If we choose 

b 

q ~ t(l-t) 2 



10 



then the tangent of (|3.4p at [z : : 0] coincides with one of the tangents 
(|3.2p . Then we have higher order intersection at [z : : 0], so that we have 
to count this node three times. Then we already have 7 intersection points. 
The remaining intersection point is a point on the surface that is in one- 
to-one correspondence with the parameter p and this gives us the desired 
parametrization. 

Taking v = t (^p + xrj J as a new parameter, we find after simple calcu- 
lation that 

bv 2 - v + a 2 b 

Z = n TT772 ( 3 - 5 ) 



(l-t)(a 2 -v 2 ) 



t(l-t) 



(bv 2 -v + a 2 b)((l - t)v 2 - 2tbv + t - (1 - t)a 2 ) 
(2bv - l)(v 2 - a 2 ) ' 



(3.6) 



The equations (|3.5p and ()3.6p parametrize the surface. They give a bijection 
with the Riemann sphere in the v- variable and we conclude once again that 
the surface has genus zero. 



3.3 The branch points 



Now we start to view (|2.6p as a four-sheeted branched covering of the Rie- 
mann sphere. That is, for each z the equation (12. 6p has four solutions for £ 
where as always we count according to multiplicity. So each z gives rise to 
four points or less on the Riemann surface. The branch points correspond to 
values of z for which there are at most three points on the Riemann surface. 
The critical t-values 

_ 1 + 2a 2 - y/l - 4a 2 b 2 _ 1 + 2a 2 + \fl - 4a 2 b 2 

^ " 2 (l + a 2 + 6 2 ) ' tc ' 2 " 2(l + a 2 + 6 2 ) (3 - ?) 

are such that 

a 2 {l-tf + b 2 t 2 = t(l-t). (3.8) 

In that case, when we put z = in (|2.6p we obtain £ 4 = 0, which means 
that the two points on the Riemann surface that correspond to z = £ = 
are both branch points. This is the critical case and we will not consider it 
any further. It is easy to see that for ab < 1/2, we have < t c \ < t c ^ < 1. 

The branch points can be found by calculating the zeros of the discrim- 
inant of (|2.6|) . With the aid of Maple, we get the following equation for the 
discriminant: 

(a 3 z 6 + a 2 z 4 + a lZ 2 + a )z 2 = 0, (3.9) 
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where 

a 3 =16aV, (3.10) 
a 2 =(-48a 2 6 2 (a - bf + a 4 + b A )t 2 

+ (96a 2 (a 2 - b 2 ) - 8a 2 )t - 48a 4 6 2 - 4a 2 , (3.11) 
ai = (48a 2 6 2 (a 4 + 6 4 - 7a 2 o 2 - a 2 - 6 2 ) 

+104a 2 6 2 - 8a 4 - 86 4 + 20a 2 + 206 2 + l) t A 

+ (48aV(4a 4 + 14a 2 6 2 + 3a 2 + b 2 ) 

+32a 4 - 208a 2 b 2 - 60a 2 - 206 2 - 2) t 3 

+ (48a 2 6 2 (6a 4 - 7a 2 6 2 - 3a 2 ) + 104a 2 6 2 - 48a 4 + 60a 2 + l)t 2 

+ (48a 4 6 2 (-4a 2 + 1) + 32a 4 - 20a 2 )t + 48a 6 6 2 - 8a 4 , (3.12) 

a =4(1 - 4a 2 6 2 )(a 2 (l - t) 2 + b 2 t 2 - t(l - t)f. (3.13) 

Note that z = is always a zero of the discriminant (|3.9p , but since z = £ = 
is a node, this is not a branch point in case t ^ i Cj i and t ^ t c ^- So (|3.9[) 
leads to the sixth degree equation 

a-iZ G + a 2 z A + aiz 2 + a = 0, (3.14) 

which has six roots. 

We show the following: 

Lemma 3.2 For every t E (0, 1), the equation (|3. 14[) /ias onZy rea/ or purely 
imaginary roots. 

(a) For t G (0, t Ci i) U (t c ,2 ; 0), there are four real roots ±z\ and ±z 2 with 
z i > z 2 > and two imaginary roots ±iz3 with z% > 0. 

(b) For t £ (t c ,i,t Ct 2), there are two real roots ±z\ with z\ > and four 
purely imaginary roots ±iz 2 and ±iZ3 with Z3 > z 2 > 0. 

Proof. Because (|3.14p only has even powers of z, it is enough to prove 
that the polynomial 

Pi(x) = a^x 3 + a 2 x 2 + a\x + oq (3.15) 

has three real roots with one negative root and two positive roots in case 
(a) and one positive root and two negative roots in case (b). 

We first examine when (|3.15|) has a multiple root. We find those values 
by looking at the zeros of the discriminant of (|3.15p . which is a polynomial 
of degree 12 in i, given by 

16i(l - t)((a + b)t - a) 2 ((a - b)t - a) 2 

x ((432a 4 6 4 + 8a 2 - 72a 2 6 2 - 1 + 8b 2 )t 2 (3.16) 
+((-432a 4 6 4 + 72a 2 6 2 + l)t(l - t) + 8a 2 (l - t) 2 + 86 2 t 2 ) 3 , 
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which was calculated with the use of Maple. We find multiple roots for 
t\ = 0, *2 = 1) ^3 = ^p5) ^4 = ^§5) an d for the zeros £5 and t& of the 
quadratic polynomial 

q(t) = (-432a 4 6 4 + 72aV + l)t(l - t) + 8a 2 (l - t) 2 + 8b 2 t 2 . 

Now it is easy to verify that 

o„2 1- * .pj.2 t 



t l-t 

takes its minimum for t E (0, 1) in t = £3 = and the minimum is 16ab. 
Thus for t G (0,1), 

q(t) > t(l - t)(-432a 4 6 4 + 72a 2 6 2 + 1 + 16a6) 
= t(l - t)(l - 2o6)(l + 6a6) 3 > 0, 

since a& < 1/2. So is and do not belong to the interval (0, 1) and the only 
value of t € (0, 1) for which (13. 15ft has a multiple root is ts = ^5. For £3 we 
find that p\ has the following roots: 

4ab(2ab+l) , (2a6 - l) 2 

xi = — 7 —75 — and £2 = ^3 ~~ 



(a + 6) 2 z 4(a + 6) 2 ' 

So for t = t% we have that pi has a double negative root and one positive 
root. 

Next, we investigate the critical i-values t = t c> i and t = t Ct 2, see (I3.7D . 
which are exactly those t 6 (0, 1) for which ao = 0. So for these values, one 
root of (|3.15p is at x = 0. We determine the sign of the other two roots by 
looking at a\. We are going to take t = t c ,i- The proof for t = t C) 2 is similar. 
If we substitute t = t c< i into (|3.12|) . we arrive at the expression 



ai = - ll^+b^Tir ( fl4((462 + + - 2a2b2 



+6 4 ((4a 2 + 1)(1 - \A " 4a 2 6 2 ) - 2aV) + 12a 4 6 4 ' 

(3.17) 

Since < ab < 1/2, it is easy to check that 2a 2 6 2 < 1 - y/l — Aa 2 b 2 . Then 
it readily follows that a\ < because of ()3.17p . Thus for t = t c> \ the 
polynomial p\ has one negative root, one positive root, and one root at 0. 
The same holds for t = t c ^- 

Now we can finish the proof of the lemma by a continuity argument. For 
t = t C) \ and for t = t Cj 2 we have three distinct real roots. Since there are 
only multiple roots for t = a /{a + b) we find three three distinct real roots 
for every t in the two intervals (0, a /{a + b)) and ifl/ifl + b), 1). Note that 
i C) i < a/ (a + b) < t c2 - 
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Since we have a double negative root and one positive root for t = 
a/ (a + b), and is only a root for t c> \ and t c ^, it follows that p± has two 
negative roots and one positive root for every t S (t Cj i,t c ^)- 

For t = 0, we find pi(x) = 4a 2 (x — a 2 ) 2 (46 2 x — 4a 2 6 2 + 1) with roots 
x\ = X2 = a 2 > and x 3 = 4:0,2 < 0. Then by continuity we find that 
there are two positive roots and one negative root for every t S (0, t c {). We 
find the same thing for t S (t c 2, 1) by looking at t = 1. 

This completes the proof of the lemma. □ 

The six roots of (|3. 14[) are all simple branch points of the Riemann 
surface. This follows from the Riemann-Hurwitz formula [241 Theorem 4.16], 
and the fact that the genus is zero. 



3.4 Sheet structure 

There are four inverse functions of (|2.6p . which behave near infinity as 



t(l- 


t) 


z 




t(l- 


t) 


b 




1-t 


+ 


b 





6M = ^-?-s + (?)' 



2z 

Uz) = --— + ^ + o(\). (3.21) 



&(*) = 7^ + ^ + 0(4), ( 3 - 2 °) 



1-t 2z 



The sheet structure of the Riemann surface is determined by the way 
we choose the analytical continuations of the ^-functions. Since all branch 
points are on the real or imaginary axis, we find that all ^ 's have an analytic 
extension to C \ C where C is the cross 

C = [-z\,zx\ U [-iz 3 ,iz 3 \. 

The analytic continuation is also denoted by £j. We consider £j on the jth 
sheet. 

Now we use the rational parametrization (|3.5|) - (|3.6|) of the surface to 
compute the corresponding curves in the u-plane. We denote the u-value 
corresponding to z £ C \ C on the jth sheet by Vj(z) and we let Uj be the 
image set of Vj(C\C). The next two figures show the regions Uj, j = 1, 2, 3, 4, 
for the cases t < i Cj i and t Cj i < t < t Cj 2, respectively. Denoting the point at 
infinity on the jth sheet by 00,-, we can determine the corresponding v- values 
from the formulas (|3.5p . (|3.6p . and (|3.18j) — (|3.21j) . Indeed we have 

z = 00 1 < — ► v = a, 
z = oo 2 « v = -a, 
z = oo 3 < — > v = 25, 

2 = OO4 < > V = OO. 
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1.5r 



1 - 




-1 - 



-1.5 1 1 1 1 ' 1 1 ' 1 

-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Figure 4: The image of the cross C = [—21,21] U [— 123,423] in the u-plane 
for < t < t Cj i. The values v = —a, v = a, and v = A are indicated with 
*. These f-values correspond to 00 on the second, first and third sheets, 
respectively. The precise values used for the figure are a = b = 0.6 and 
t = 0.25. 



These four u-values are independent of t, and they are indicated with * in 
Figures [U and [5] (except for v = 00, which is of course not visible). Note 
that —a < a < ^. 

The figures show what the analytic continuations are of the functions £j 
and how the sheets of the Riemann surface are connected. First of all, we 
note that any part of the boundary of ujj that is not a part of the boundary 
of some u>k with k 7^ j, corresponds in the z-plane to a part of the cross C 
where £j has equal boundary values from both sides. This part of the cross 
is then removed from the cut on the j'th sheet. 

In this way, we see that for t < t Cj i, £1 is defined with a cut [22, 21] only, 
£2 with a cut [—21,-22], £3 with two cuts [22,21] and [—^23,^23] and £4 is 
defined with two cuts [—22,-21] and [—223,423]. The sheets are connected 
as shown in Figure [H 

For t c 1 < t < t Cj 2 we have the two pairs of purely imaginary branch 
points ±422 and ±423, where we choose 23 > 22. In the situation of Figure 
we have that £1(422) = £2(^22) and £3(423) = £4(423). This is always the 
case if i Cj i < t < a/ (a + b). Then we have that £1 is analytic with a cut 
on [0,2i] U [—422,422], £2 has a cut on [— 21, 0] U [—422,422], £3 has a cut 
on [0,2i] U [—423,423], and £4 has a cut on [— 21, 0] U [—423,423]. The sheet 
structure is then as shown in Figure [71 

For t = a/(a + b) we have that 22 = 23 and for a/ (a + b) < t < t c ^ 



15 



Figure 5: The image of the cross C = [—zi,z{\ U [—iz^^iz^] in the u-plane 
for i C)1 < t < t Cj 2- The values v = —a, v = a, and v = ■% are indicated with 
*. These v-values correspond to oo on the second, first and third sheets, 
respectively. The precise values used for the figure are a = b = 0.6 and 



we find that the role of Z2 and Z3 are reversed, so that then the first and 
second sheets are connected along [— iz^jiz^] and the third and fourth along 
[-iz 2 ,iz 2 \. 

For t C) 2 < t < 1, we have a similar sheet structure as in Figure O except 
that the cut on the vertical segment [— izs,izs] now connects the first and 
the second sheets. 

3.5 Properties of ^ 

The sheet structure induces jump relations between the ^-functions along 
the cuts that we will use later on. At the branch point Z\ we have for a real 
constant c\ > that 



In case t £ (0,i Cj i) U (t C} 2, 1) we also have that at the branch point z 2 there 
is a real constant c 2 > such that 



t = 0.45. 



£1(2) = £1(21) + c i( z ~ z\) 
£30) = £l(zi) - ciO - zi) 



/2 + 0(z _ Zl)) 

l 2 + 0(z- Zl ). 



(3.23) 
(3.24) 



£1(20 = £1(^2) - c 2 (z 2 - z) 

£3(*0 = £1(^2) +C 2 (^2 - 2) 



l/2 +0 ( Z2 _ z)) 
V2 +0(Z2 _ Z) . 



(3.25) 
(3.26) 
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Figure 6: The sheet structure of the Riemann surface for ab < 1/2 and 
< t < t cA . 



Here the square root is taken with a branch cut along the negative real 
axis. Similar behavior also holds for the functions £2 and £4 near the branch 
points —Z\ and — Zi. 

From the sheet structure we also see how the ^-functions are continued 
along the cuts. We will not list those relations explicitly, but they will be 
(tacitly) used in the future. 

What is also important is that for any closed contour 7 that does not 
intersect the cut on the j'th sheet, we have 



This follows from the fact that the residues at infinity for each of the £j 
functions is ±±, see ([535]> - (|5^I]) . 




(3.27) 
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Figure 7: The sheet structure of the Riemann surface for ab < 1/2 and 

t £ (t c ,l,t Cj 2)- 

3.6 The Xj functions 

The A- functions are primitive functions of the £- functions. We define them 

as 



Ai 


» = 


/ £i(s)ete, 


A 2 


>) = 


/ 6(s)ds + c, 

7-21 + 


A 3 


>) = 


/ C3(s)ds, 


A 4 


» = 


I ii{s)ds + c, 

J—Zl — 



(3.28) 



where c is chosen so that Xs(iz3) = X^izs) (if < t < a/ {a + 6)), or 
Ai(i^3) = X2(iz%) (if a/(a + b) < t < 1). The path of integration in each 
integral in (^28l) is in C \ (C U (-oo, -z{\). 

From these definitions, and from the fact (|3.27p that the periods of the 
^-functions are half integers, it then follows that e nX is analytic on the 
Riemann surface (recall n is even). That is, if Cjk is the cut connecting the 
jth. and kth sheets, then 
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The behavior as z — ► oo follows from (|3.28|> and (|3.18p - (l3.2ip and is 
given by 



z 2 az 1 , , „ / 1 



Al(3 » = 2«(r^)-T-2 ta2+!l + UJ' (3 - 30) 

. / s -z 2 1 , „ /1\ . 

^ = iR + i-2 hz+i2+o y- (3 - 3i) 

\3(z) = -^- + llnz + k + o(-) , (3.32) 



1-f 2 V^, 

X ^ = ~Y~t + \ lnz + k + °(S) ' (3 ' 33) 

where h,h,h^4 are certain integration constants. 

4 Steepest descent analysis for case < t < t c i 

We will do the steepest descent analysis in some detail for the case < t < 
t Ct i. The case t c> 2 < t < 1 follows from this by symmetry of the problem. 
For the case i Cj i < t < t Cj 2, we refer to Section O 

The steepest descent analysis is based on the one given in [JJ. A major 
role is played by the Riemann surface, which for for < t < t C) i, has the 
sheet structure as shown in Figured! 

4.1 First transformation: Y \— > U 

In the first transformation we normalize the RH problem at infinity. We 
define 



(4.1) 



where L is the constant diagonal matrix 

L = diag (V* 1 , e~ h , e~ h , e"' 4 ) . (4.2) 

Here the constants L, j = 1, 2, 3, 4, are the constants that appear in (j3.30jl - 
([3331 - Now U is defined and analytic in C \ (R U [-iz 3 ,iz 3 ]). Then ([2TT4]) , 
(|3.30p - ()3.33p . and (|4.1|) imply that U is normalized at infinity: 

U(z) = l + 0(l/z) asra^oo. (4.3) 
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Using (14, ip and (12.130 we find the jumps for U. On the real line, we get 

e n(A 3+ -Ai_) e n(A 4 +-Ai_)\ 

g n(A4+— A 2 _ 
3 n(A 3+ -A 3 _) q 



U+ = U- 



/ e n(Ai+-Ai_) 







e n(A 2 +-A 2 _) e n(A 3+ -A 2 _) e n(A 4+ -A 2 _) 







n(A 4+ -A 4 _) J 



(4.4) 



while on the vertical segment [—iz3,izs], which is oriented upwards (so that 
U+{z) (U-(z)) is the limiting value of U(z') as z' — > z from the left (right) 
half-plane) we have the jump 



U+ = C/_diafffl,l,e n(A3 +- A3 - ) ,e n(A4 +- A4 - )N l . 



(4.5) 



Because of (|3.29j) we can simplify the jump (|4.4p on the various parts of the 
real line. The result is the following RH problem for U : 

RH problem for U 

(1) U is analytic onC\(RU [— 123, 223]). 

(2) U satisfies the following jumps on the real line: 



U 4 



f\ e "(A 3 +-Ai_) e n(A 4+ -Ai_)\ 

1 e rl (- > *3+-A 2 _) g n(A 4+ -A 2 _) 

1 

\0 1 / 



on (—00, —z\), 
(4.6) 



U+ = U- 



(\ e n(A 3+ -Ai_) e n(A 4+ -Ai_)\ 

e «(A 2 +-A 2 _) e n(A 3+ -A 2 _) y 

1 

\0 e n(A 4+ -A 4 _)y 



U+ = U- 



/l e n ( A 3+-Ai_) e n(A 4+ -Ai_)\ 
1 g n.(A 3 +-A 2 ) e n(A 4+ -A 2 ) 





\o 



1 






1 



/ 



U+ = U- 



/ e n(Ai + -Ai_) q 1 e n(A 4 -A!_)^ 

1 e ™(^3+-A 2 ) e "(A 4 -A 2 ) 

e «( A 3+- A 3-) 

1 / 



on (-zi, - 
(4.7) 

on (-z 2 ,0)U(0,z 2 ), 
(4.8) 

on (z 2 ,z 1 ), 
(4.9) 
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U+ = U- 



(\ e n ( A3 ~ Al ) 

1 

\o 



e ra(A4— Ai)\ 
e n(A 4 -A 2 ) 



1 



on 0i,+oo). (4.10) 



Finally, on the vertical segment, which is oriented upwards, 

U+ = L/_diag ^l ) l )e «(^ + -A3-) e n(A4+-A4-)^ Qn (_ iza izs y 

(4.11) 

(3) U(z) = 1 + 0(l/z) as z -> oo. 

Now it turns out that not all entries in the jump matrices for U are well- 
behaved as n — > oo. Ideally, one would like to have exponentially decaying 
off-diagonal terms e n ^ j+ ~^ k - s> in all of the jump matrices in (I4,6p - (|4.1ip . 
and oscillating diagonal entries e n ^ j+ ~ Xj -" > . However, this is not the case in 
the present situation. For example, the entry e n ( A3 + -A3_ ) in the jump matrix 
in (|4.1ip is exponentially increasing as n — ► oo, since Re A 3+ > Re A3_ on 
the vertical segment (— 1Z3, izs). 



4.2 Second transformation: U i— > T 

We use a trick to remove the exponentially increasing entries. This involves 
the global opening of a lens, which was introduced in [3], and also used in 

We take a closed curve S, consisting of a part in the right half-plane 
that connects — iz?, and iz$, that is symmetric with respect to the real axis, 
and that satisfies 

Re X 4 (z) < Re X 3 (z) for z £ E n {Re z > 0}. (4.12) 

The curve E intersects the positive real line in a point x* and we assume 
that x* is sufficiently large so that 

ReA 4 (x) < ReA 3 (x) < ReAi(x) < ReA 2 (a;) for x > x* , x G R. 

The part of E in the left half-plane is the mirror image with respect to the 
imaginary axis. Then 

Re A 3 (z) < Re \ A {z) for z G E n {Re z < 0}, (4.13) 

and 

ReA 3 (x) < ReA 4 (x) < ReA 2 (x) < ReAi(a;) for x < -x*, x G R. 

Note that the behavior (|4.12|) - (|4.13|) is valid near infinity because of (|3.32|) - 
(|3.33p and we have that the inequalities remain valid in a domain whose 
boundary contains the branch points Hz^. 
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Figure 8: The behavior of the real part of Xj — Xk when a = b = 0.6 and 
t = 0.05. Only some relevant curves are shown where two real parts coincide. 



This is illustrated in Figures [8] and [9] where the solid curves are the 
curves where Re A3 = ReA4. We have three such curves that emanate at 
equal angles from IZ3. One of them extends to infinity along the positive 
imaginary axis. The other two continue to the branch point — iz^ and to- 
gether form a closed contour. This closed contour can either contain the 
intervals [—Z\, —Z2] and [z2,Zx] in its interior as in Figured! or it intersects 
these two intervals as in Figure [HJ 

For later convenience, we also choose £ as the analytic continuation of 
the closed contour where Re A3 = ReA4 in a neighborhood of Hz^. For an 
example how to choose S, see Figure [TU1 

Now we are ready to introduce the second transformation, see also [U 
section 4]. We define T as 



T = U 



T = U 



outside E, 

/l 

1 

1 

\0 - e n ( A 3- A4 ) 






1/ 



(4.14) 



inside E in the left half-plane, 



(4.15) 
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u 



/I 











\ 





1 
















1 


_ g n(A 4 - 


-As) 


\o 








1 


/ 



inside E in the right half-plane. 



(4.16) 

Then T has jumps on the contours shown in Figure [IT] and a straight- 
forward calculation shows that T satisfies the following RH problem: 

RH problem for T 

(1) T is analytic on C \ (R U [— iz%, iz%] U E). 

(2) On the real line, T has the following jumps: 

/l e n ( A3 +~ Al ~) e n(A 4 +-Ai_)\^ 



1 

\0 



s n(A 3+ -A2-) n(A 4+ -A 2 _) 



1 





/ 



(I e «(A4+-Ai_)\ 



1 
1 
\0 



e "(A4+-A 2 -) 



1 



on (— oo, —x*), 

(4.17) 

on (-x*,-zi), (4.18) 
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-0.8 1 1 ' 1 1 1 1 

■1.5 -1 -0.5 0.5 1 1.5 

Ftez 

Figure 10: Opening of the global lens £ that connects ±iz$. It is such that 
Re A3 > Re A4 on X in the right half-plane, and Re A4 > Re A3 on £ in the 
left half-plane. Here we have taken a = b = 0.6 and t = 0.05. 



T+ = T_ 



T + = T_ 



/10 

e n(A 2+ -A 2 -) q 

1 










Vo 



e n(A4 + -Ai_)\ 



1 



e n(A4+— A4_ 



V 



/1 

1 



\o 

/l 

1 



\o 

/ e n(Ai+-Ai_) 




V 



g n(A 4+ -Ai_)\ 



on (-zi, -z 2 ), 

(4.19) 



e n(A4+— A2) 



1 

e n(A 3+ -Ai_) 
e n(A3+— A2) 

1 



on(-z 2 ,0), (4.20) 



/ 






1/ 



on 



(0,^2), 



(4.21) 



1 



e re(A 3+ -A 2 ) 
e n(A3 + -A3_) 








1/ 



on (z 2 ,^i), 



(4.22) 
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Figure 11: Jump contour for the RH problem for T. The contour E has 
clockwise orientation. The 4x4 matrix valued function T is analytic outside 
these contours, and has jumps (|4.17|) ~ (|4.27|) along the various parts of the 
contour. 



/I 



1 e n ( A3 ~ A2 ) 

1 

\0 1/ 

(I e "( A 3-^i) 

1 g"(A3-A2) 

1 

\0 



on (zi, x 



(4.23) 



e n(A4-Ai)\ 
e n(A4-A 2 ) 



1 



on (x*,+oo). (4.24) 



On the vertical segment (with upwards orientation), we have 

\ 





1 

= n(A 4+ -A 4 _) 



on (-iz 3 ,iz 3 ). 



(4.25) 



/l 

10 



\0 -1 e n(A4 +~ A4 -V 

The jumps on E are as follows, where we take clockwise orientation 
on E, so that T + {z) (T_(z)) for z £ S, is the limiting value of T(z') 
as z' — » z from outside (inside) of E: 

0\ 



/I 

1 



\0 

/l 

1 



e n(A 3 -A 4 ) 

\ 



1 e n ( A4 ~ A3 ) 
\0 1 ) 



on {z G E | Rez < 0}, (4.26) 



on {z £ E | Rez > 0} (4.27) 
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Figure 12: Opening of the lenses T\ and T2 around the intervals [—zi, —Z2] 
and [z 2 ,zi\. 



(3) T(z) = 1 + 0(l/z) asz^oo. 

Now it follows that the jump matrix in (14.250 on the vertical segment 
( — iz$, 1Z3) tends to a constant matrix since Re A4+ < Re A4_. This may also 
be deduced from Figures [8] and EE since we have Re A3 < Re A4 in the right 
half-plane within the closed contour in these figures. The strict inequality 
remains valid up to the vertical segment (— izs, IZ3) so that Re A3_ < Re A4_, 
since the orientation of the segment is upwards. Since Re A4± = Re A3zp the 
inequality ReA4 + < ReA4_ follows. 

All the jump matrices on the real line tend to the identity matrix as 
n — > 00, except for the ones on the intervals (— Zi,— z%) and (22,21). The 
1,1-entry and the 3,3-entry of the jump matrix on (22,21), see (14.221) . are 
rapidly oscillating for large n, and the same holds for the 2,2-entry and the 
4,4-entry of the jump matrix on (— £1,— £2), see (14.190 . These oscillating 
entries are turned into exponentially decaying ones in a standard way by 
opening lenses around the intervals [— £1,— £2] an d [22,21]. This will be 
done next. 



4.3 Third transformation: T 1— > S 



We now turn the oscillatory entries on the diagonal of (14.190 and (|4.22p 
into exponentially decaying ones. This is done in a standard way by open- 
ing lenses T\ and T2 in (small) neighborhoods of (— z%,— 22) an d (22,21), 
respectively, as shown in Figure [121 Define S as follows: 



S = T outside Ti and T2, 



/I 








1 





A 4 ) 






1/ 



(4.28) 



in the upper lens region around [—Zi, —Z2] 



(4.29) 
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Figure 13: Jump contour for the RH problem for S. The contour £ has 
clockwise orientation and the upper and lower lips of the lenses are oriented 
from left to right. The 4x4 matrix valued function S is analytic outside 
these contours, and has jumps (|4.33p - (|4.45p along the various parts of the 
contour. 



S = T 



S = T 



S = T 



(\ 0\ 

1 

10 

\0 e < X2 ~ x ^ l/ 



( 1 0\ 

10 

_ e n(Ai-A 3 ) o 1 o 

V o 001/ 



/ 1 0\ 

10 

e n(Ai-A 3 ) 10 

V o ooi/ 



in the lower lens region around [— z±, —Z2], 

(4.30) 

in the upper lens region around [z 2 , zi], 

(4.31) 

in the lower lens region around [Z2 , z\\. 

(4.32) 



Then S satisfies the following RH problem on the contour shown in 
Figure 13 

RH problem for S 

(1) S is analytic on C \ (R U S U T x U T 2 U [-iz 3 , iz 3 }), 
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(2) On the real line, S has the following jumps: 



/l e n ( A 3+~ A i-) e n(A 4 +-Ai_)\ 

1 



Vo 



S + = S- 



S + = S- 



(I 






S+ = S- 







/l e «(A4+-Ai- 



1 
1 
\0 



on (— 



e n(A4 + -A2-) 



1 / 

e «(A4+— Ai 




e n(A2+— Ai_) 



■v*,-zi), 



on (— oo, —x*), 
(4.33) 

(4.34) 








Vo 

/l 

1 

1 

\0 



1 





e n(A4 + -Ai_)\ 
e ra(A 4+ -A 2 ) 



1 J 



1 






on (-zi, -z 2 ), 
(4.35) 

on(-z 2 ,0), (4.36) 



S+ = S- 



S+ = S- 



S + = S- 



S+ = s. 



/I 







1 





e n(A 3+ -Ai_) 
e n(A 3+ -A 2 ) 

1 



\ 






on (0,2:2), 



(4.37) 



Vo 







1/ 






/ 







1 




o\ 




e n(Ai + -A 2 ) 1 


e ™(A3+ 


-A 2 ) 







-1 












V 












1/ 


( l 


e "(A3-Ai) 


o\ 











I e n(A 3 -A 2 ) 














1 







on (zi, x 


\o 





1/ 








/I 


e «(A3-Ai) 


e n 


(A 4 - 


Ai)\ 







I e "(A 3 -A 2 ) 


e™ 


(A4- 


A 2 ) 







1 









c 


\o 







1 


y 





on (z 2 ,z 1 ), 



(4.38) 



(4.39) 



on(x*,+oo). (4.40) 
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On the vertical segment, we have that 



(1 \ 

10 

1 

\0 -1 e n(A 4+ -A 4 -)y 



on (-iz 3 ,iz 3 ). 



(4.41) 



On E, the jumps of S are 



5+ = S. 







/l 

1 

1 

\0 e n ( A3 " A4 ) iy 



fl \ 

10 

1 e "( A 4- A 3) 

\0 1 y 



on E in the left half-plane, 

(4.42) 

on E in the right half-plane. 

(4.43) 



The jumps on the lenses r*i and T2 around [—zi, —Z2} and [Z2, z\] are 

(4.44) 



S-i 



/l 





0\ 

1 
1 

\Q e n(A 2 -A 4 ) l J 

(I 0\ 

10 

e n(Ai-A 3 ) 10 

V 001/ 



on Ti 



on To 



(4.45) 



(3) S = I + 0(l/z) asz^oo. 

Now it can be checked that all non-constant entries in the jump matrices 
for S tend to as n — > 00. This could be done by completing the Figures [8] 
and [9] by drawing all curves where Re Xj = Re for j ^ k from which one 
can deduce the sign of Re (Xj — A&) in various regions in the complex plane. 

For the jumps on (— z±,— z?) and (z2,z±) one can also use a Cauchy- 
Riemann type argument. For the interval (z2,z\) this is based on the fact 
that (Ax — A3)+ is purely imaginary on (z2,z\) and its derivative is (£1 — 
£3)+ = 2Tm£i + where Im£i + > 0. Then by the Cauchy-Riemann equations 
we have that Re (Ai — A3) < in a neighborhood above (Z2, z±). 
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4.4 Parametrix away from the branch points 

We are now going to solve the model RH problem, where we ignore all 
exponentially small entries in the jump matrices for S. So we only keep the 
jumps on (—21,-^2); ( z 2,zi), and (— 1Z3, IZ3) (which are also the cuts for 
the Riemann surface) and we look for a 4 x 4 matrix valued function iV that 
satisfies the following RH problem: 

RH problem for iV 

(1) N is analytic on C \ ([—21,-22] U [^2,^1] U [ — ^3,^3]), 

(2) N satisfies the following jumps along the cuts: 



N+ = iV_ 



iV_ 
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iV_ 
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on (-zi, -z 2 ) 



on (z 2 ,zx) 



on 



(4.46) 



(4.47) 



(4.48) 



(3) as z — > 00, we have that 
iV(z) — J + 0(l/z). 



(4.49) 



By means of the rational parametrization of the Riemann surface we can 
transform this problem to the complex f-plane and thereby solve this RH 
problem explicitly. 

Recall that we have the mapping v = Vk(z) that maps the kth sheet of 
the Riemann surface to the part uik of the f-plane, see also Figure [H We 
look for a solution for the RH problem for iV in the form 



N{z) = {F 3 {v k {z))) 



j,k=l,...,A 



where Fj, j = 1, ... ,4 are four functions on the f-plane. Then the jump 
conditions on iV are satisfied if the Fj are analytic except for the parts of the 
boundaries of the domains lo^ that are indicated in Figure [T4"l and on these 
parts we have Fj + = —Fj-. The asymptotic condition on N is satisfied 
provided that 



Fj(v k (oo)) = S jk , 



for j,k= 1,. . . ,4. 
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Figure 14: Jump contours in the v-plane for the scalar RH problems for the 
functions Fj, j = 1, 2, 3, 4. 



For each j, we then have a scalar RH problem for Fj that can be solved 
explicitly with elementary functions. 

The precise form of Fj is not important for what follows, but we do need 
that Fj(v) behaves like (v — vq)~ 1 / 2 as v — » v o and vq is one of the endpoints 
of the jump contour in Figure [TH For N this implies that 



where zq is any of the branch points izi, ±£2, ±^3- 
4.5 Parametrix near the branch points 

The jump matrices on N and S are not uniformly close to each other near 
the branch points. That is why we need to treat these points separately, 
and construct a local parametrix P around these branch points. 

We are going to construct a local parametrix around Z\ . The local para- 
metrices around —Zi, ±2*2, ±123 can be found in a similar way, and are 
therefore not further discussed here. Consider a small but fixed disc Ug 
with radius d around Z\ that does not contain z%. We then look for a 4 x 4 
matrix valued function P such that 

RH problem for P around z\ 

(1) P is analytic for z E Ug \ (T2 U K), for some 6q > 5, 




as z 
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(2) P has the following jumps on the real line: 



P+ = P- 



P+ = P- 



( 1 0\ 
10 
-10 

V o 001/ 

/l e n ( A3 ~ Al ) 0\ 

1 

1 

\0 1/ 



on U s n [z 2 ,zi], 



on Us n [z%, +oo) 



(4.50) 



(4.51) 



and on the lens T 2 around [z 2 , z\], we have 



P+ = P- 



( 1 0\ 

10 

e n(Ai-A 3 ) 10 

V 1/ 



on u s nr 2 . 



(4.52) 



(3) As n — > oo, 

P(z) = iV(^)(J + 0(l/n)) uniformly for z £ dU 5 \ 



ur 2 ). 

(4.53) 



Note that the jumps for P are not exactly the same as the jumps for S. 
They differ by the entries e n ( A3 +~ A2 ) and e n ( Al +~ A2 ) which are exponentially 
small as n — > oo. 

The RH problem is solved in a standard way with Airy functions |15l 
[161 HE]. We have that 



i (Ai-A 3 )(2) 



2/3 



(4.54) 



is a conformal map that maps a neighborhood of Z\ onto a neighborhood of 
the origin such that fi(z) is real and positive for z > z\. We open the lens 
around [z 2 , Z\\ such that f\ maps the part of T 2 in this neighborhood of z\ 
to the rays with angles ^ and — respectively. We put 

y (s) = Ai( S ), yi ( S ) = e 2 ^ 3 Ai(e 2 -/ 3 S ), y 2 ( S ) = e" 2 ^ 3 Ai( e ~ 2 -/ 3 S ), 

(4.55) 
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where Ai is the usual Airy function. Define the matrix ^ by 



*00 



/yo(s) -y 2 (s) 0\ 

10 

yo( s ) o -y 2 (s) 

V o o o 1/ 

-y 2 (s) 0\ 
10 
-y' 2 (s) 

V o o o 1/ 

/-y 2 (s) o yi (s) o\ 

10 

V o 001/ 

/y (s) j/!( S ) 0\ 
10 
yo(s) y[(s) 

V o o o 1/ 



for args G (0, 27r/3), 



for args G (27r/3,7r) 



(4.56) 



(4.57) 



for args e (-vr, -2vr/3), (4.58) 



for args G (— 27r/3, 0). 



(4.59) 



Then, for any analytic prefactor E, we have that 

P(z) = E(z)V (n 2 / 3 diag (e^ x ^~ x ^\ 1, e ~f (M*)-^)), l) 

(4.60) 

satisfies the parts (1) and (2) of the RH problem for P. If we choose E as 



^■-v^v ".. : " : ;: : „ 1/4 (4.ei) 

o 1 1/ 

then £" is analytic and the part (3) is satisfied as well. 
4.6 Fourth transformation: S i— > R 

In the final transformation, we define the matrix valued function R as 

R(z) = S(z)P(z)^ 1 in (small) discs around ±z±, ±z 2 and Hz^, 

(4.62) 



/I 





-1 0\ 










o\ 





1 








1 








— i 





-i 








n -vvr 1/4 





\o 





1/ 


V o 








V 



i?(z) = S(z)N(z) 



-i 



outside the discs. 



(4.63) 



Then 72 is defined and is analytic (more precisely, has analytic continuation 
to the region) outside the contour shown in Figure [I5j 
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Figure 15: Jump contour for the RH problem for R. The 4x4 matrix 
valued function R is analytic outside these contours, and has jumps R + = 
R_(I + 0(l/n)) uniformly on all parts of the contour. 



From the matching condition (|4.53|) (and similar ones around the other 
branch points), it follows that on the circles around the branch points there 
is a jump 

R + = R-(I + 0(l/n)) uniformly as n — > oo. (4.64) 
On the remaining contours, the jump is given by 

R+ = R-(I + 0(e~ cn )) asn^oo (4.65) 
for some constant c > 0. Together with the asymptotic condition 

R(z) = I + 0(l/z) asz^oo, (4.66) 
it then follows as in [2Q E5 ESI ESI [22] that 

R( z )=I + o(-— r -) asn^oo, (4.67) 

\n{\z\ + l)J 

uniformly for z £ C\ Tr, where Tr is the jump contour for the RH problem 
for R, see Figure EH 

5 Proofs of theorems for case < t < t Cj \ 

Having (|4.67p we can now prove Theorems I2.1H2.3I in the same way as in [6] . 
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5.1 Proof of Theorem ED 



Take x,y £ (z2,Zi). We follow the transformations Y i— > U i— > T ^ S, to 
obtain from (|2.16p that 

n(h(y)-h(x)) 

K n (x,y) = - |_ginImAi + (s() q e -mImAi + (y) e n(A 4 (j/)-Re Xi + (y))\ 



2iri(x — y) 
xS^(y)S + (x) 



/ e -inlvn X 1+ (x)\ 




\ 







(5.1) 



/ 



where 



M^O = o ( ReA i+( x ) + ReA 3+(x)) 



for x > 0. 



2 v—-—— ■ -----oi-w/ 2(1 -i) 

As in pU Section 9] we obtain from (|4.67|) that 

S^ 1 (y)S + (x) = 1 + 0(x - y) as y -> x 

uniformly in n, and therefore we get that 

^ ) = e n(h(y)-h( x) ) ( sin(nlm (A 1+ (x) - Ai+(y))) + 

V vr(x-y) 

where the 0(1) holds uniformly in n. Letting y — > x we find 



K n (x,x) 



n 



7T 



lm£ 1+ (x) + 0(l). 



(5.2) 
(5.3) 
(5.4) 



p(x) 



(5.5) 



For x 6 (— Z2, —z\) we get the similar relation but with £i + (x) replaced by 
^2+(x). Thus (12. 4h holds with p defined by 

-Im£i+(x) for x > 0, 

^Im^2+(^) for x < 0. 

The further statements in Theorem 1 2 . 1 1 are now easy consequences of the 
properties of £i and £2- 

5.2 Proof of Theorem Q 

Let xo G (z2,zi) and take 

s = + — , , y = x + 



np(xo)' 7 ~ u ' np(x ) 
Then for n large enough, we have x, y £ (z2,z±), so that (|5.ip holds. Thus 



np(xo) ' ?!"(« — v) 

This leads to (pTTUj) as in [6 S section 9.2]. 
The proof for xq £ (—^2, — 21) is similar. 
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5.3 Proof of Theorem 12^1 



Take x = z\ + (ci ^ )2/3 and y = z\ + 
the transformations T ^ 



If u, u < 0, then we can follow 



(cm) 2 / 3 ' 
*S* i ^ R^ to find from (|2TT6j) 



1 



(cin) 2 / 3 



K n (x,y) 



2ixi{u — v) 



(-1 1 



nRe(A 4+ (j/)-Ai + (y)) 



) 



x [n 2 l' i f 1 {y))E- l {y)R- l (y)R + {x)E n {x) 



x^+(n 2 / 3 /i(^) 



n\ 



1 

w 



(5.6) 



The term e nRe is exponentially small, and does not contribute 

to the limit. Then we can use the arguments of [6j section 9.3] to obtain 
(|2.1ip . Similar arguments give (|2.1ip in case u > and/or v > 0. 
Likewise we get (|2,12p . 



6 Steepest descent analysis and proofs of theo- 
rems for case t Cj i < t < t Cj 2 

The steepest descent analysis is somewhat different for the case t Cj i < t < 
t c> 2, due to the different sheet structure of the Riemann surface, see Figure El 
However the main lines in the proof remain the same. We only point out 
that now we have two pairs of purely imaginary branch points, ±iz2 and 
±izs. In the second transformation we need to open two global lenses in 
order to remove the exponentially increasing entries in the jump matrices. 
After the small opening of a lens around (— zx,z±) we then construct local 
parametrices at the branch points with Airy functions. 

At the end of the transformations we arrive at a RH problem for R with 
jumps on the contour shown in Figure [16J The jump conditions for R are 
R + = i?_(7 + 0(l/n)), uniformly on all parts of the contour. Then similar 
arguments and calculations lead to the proofs of Theorems I2.1H2.3I for the 
case t c> i < t < t C) 2- We refer to [H] for complete details. 
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